Particle diode: Rectification of interacting Brownian ratchets 
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Transport of Brownian particles interacting with each other via the Morse potential is investigated 
in the presence of an ac driving force applied locally at one end of the chain. By using numerical 
simulations, we find that the system can behave as a particle diode for both overdamped and 
underdamped cases. For low frequencies, the transport from the free end to the ac acting end is 
prohibited, while the transport from the ac acting end to the free end is permitted. However, the 
polarity of the particle diode will reverse for medium frequencies. There exists an optimal value of 
the well depth of the interaction potential at which the average velocity takes its maximum. The 
average velocity v decreases monotonically with the system size N by a power law v oc iV -1 . 
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I. INTRODUCTION 

Brownian motors [l|, rectifying nonequilibrium ran- 
dom walks on asymmetric potentials, have recently been 
proposed for a variety of applications, including in bio- 
logical systems @, as well as their potential technolog- 
ical applications ranging from classical non-equilibrium 
models fS] to quantum systems Ratchets have been 
proposed to model the unidirectional motion driven by 
zero-mean non-equilibrium fluctuations. Broadly speak- 
ing, ratchet devices fall into three categories depending 
on how the applied perturbation couples to the substrate 
asymmetry: rocking ratchetsjH, flashing ratchetsQ, and 
correlation ratchets Q. Additionally, entropic ratchets, 
in which Brownian particles move in a confined structure, 
instead of a potential, were also extensively studied Q. 

In many physical situations, one deals with not a sin- 
gle Brownian particle, but rather the with arrays of a 
finite number of interacting particles, which are usually 
in a periodic potential and acted upon by some external 
force. There has been a lot of interest in recent years 
in the dynamics of interacting Brownian particles (9-fl4|. 
The reason for this interest is twofold [15j . First, ex- 
periments have provided a wealth of information about 
the motion of individual colloidal particles. A system 
of interacting Brownian particles is the simplest model 
of a colloidal suspension. Second, interacting Brown- 
ian particles constitute the simplest model system, on 
which one can test techniques and approximations of 
nonequilibrium statistical mechanics. Many-particle sys- 
tems may exhibit some features not found in the single- 
particle counterparts, such as phase transition s, sp onta- 
neous ratchet effects, and negative mobility |9l-tl4|. The 
interacting Brownian ratchets been proposed for a vari- 
ety o f ap plications, including molecular motors 0], fric- 
tion [l||, diffusion of dimers on surfaces [ll}, diffusion 
of colloidal particles [l8[ , DNA translocation through a 



nanopore [19j, charge density waves [20j], and arrays of 
Josephson junctions [2l[ . 

In the present work, we study the transport of Brow- 
nian particles interacting with each other via the Morse 
potential by applying an ac driving force at one end of 
the chain. From the numerical simulations, we find that 
the system can behave as a particle diode for both over- 
damped and underdamped cases, the transport is per- 
mitted when the ac driving force is applied at one end 
of the chain, while the transport is prohibited when the 
force acts on the other end of the chain. 



II. MODEL AND METHODS 

We consider an open system of N point-like Brown- 
ian particles interacting with each other via the Morse 
potential V, with the substrate through the asymmetric 
periodic potential U. As to obtain the directed trans- 
port, a sustained time-periodical force F(t) is applied at 
one end of the chain. The equations of motion for the 
underdamped case are 

mil =f(x l )-G(x 2 -Xx)-'yxi+€ l (t)+6 q , 1 F(t), (1) 



mx l = f(xi) + G(xi - x i - 1 )-G(x i+ i - x { ) - jx\ + & (t) , 

(2) 

mx N = f(x N ) + G(x N -x N -i)-jx N +£ N (t) + 6 q . N F(t), 

(3) 

where xi is the position of the ?th particle, m is the mass 
of the particles, and N is total number of the particles. 
f = — W- and G = — 7 is friction coefficient and the 

J ox ox 1 

noise terms satisfy the fluctuation dissipation rela- 
tions )) = 2D5(t — t). D is the noise intensity. 
The dot stands for the derivative with respect to time t. 
Where q is a nature number. For q = 1 and q — N, the 
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where k is well depth of the interaction potential and 
a controls its width. The interaction potential is asym- 
metric and the repulsive force is dominating. F(t) is a 
sustained time-periodical force, 



F(t) = A sm(ut), 



(9) 



where Aq is the amplitude of the external force and uj is 
its frequency. 

Here we mainly focus on the transport of the driven 
particles. The average velocity v j of the ith particle is, 



Vi(t) = lim 



(Xj(t) - Xj(t )) 

t-t 



(10) 



where to and t are the initial and the end times for the 
simulations, respectively. From Fig. 2 we can find that 
after the system reaches a stationary state, Vi is indepen- 
dent of site position i. Therefore, we use the asymptotic 
velocity of the particle N/2 + 1 to measure the transport 
of the system v — Ujv/2+i- 



FIG. 1: (Color online) (a) Schematic illustration of the model. 
Case A (q = 1): the driving force acting on the left end. Case 
B (q = N): the driving force acting on the right end. The on- 
site potential U(x) is asymmetric [shown in Eq. (7)]. The ac 
driving force F(t) = Ao8sn(ut). U = 1.0 and A = 1.0. Note 
that it is possible that two particles share a site or that the 
chain is torn apart, (b) The interaction potential V(x) is a 
Morse potential [shown in Eq. (8)] . The interaction potential 
is asymmetric. The parameters are k = 1.0 and a = 1.0. 



ac driving force is applied at the left end and the right 
end, respectively. 

When the inertia is negligible compared to the viscous 
damping, Eqs. (l)-(3) can be rewritten for overdamped 
case, 



7x1 = f{x x ) - G(x 2 - x x ) + &(t) + S q>1 F(t), 



(4) 



— •— m =0.0145 

— *— (o=0.145 
— A— t0=1.0 
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FIG. 2: (Color online) The dependence of average velocity 
V on the site i for overdamped case. It is found that each 
particle has the same average velocity. The other parameters 
are k = 1.0, A = 2.0, A = 1.0, and D = 0.5. 



jxi = f(xi) + G(xi - Xi_i) - G(x i+1 - x{) + £i(t), (5) 

JXn = f(x N ) + G{x N - x N -i) + ^jv(*) + S qiN F(t). (6) 

U (x) is an asymmetrically periodic potential [shown in 
Fig. 1(a)] 

U(x) = -U [mn(x) + j sm(2x)], (7) 

where Uo denotes the height of the potential and A is its 
asymmetry parameter. 

The interaction potential between the nearest-neighbor 
particles is given by the Morse potential [shown in Fig. 
1(b)], 

V(x) = k[l -exp(-ax)] 2 , (8) 



In our simulations, the equations of motion are inte- 
grated by using the second order Stochastic Runge-Kutta 
algorithrn|22l] with the small time step At = 10 ~ 3 . All 
quantities of interest are averaged over 500 different re- 
alizations. The simulations are performed long enough 
to allow the system to reach a nonequilibrium steady 
state in which the local velocity is a constant along the 
chain. To obtain a steady state, the total integration is 
typically 10 7 time units. We have checked that this is 
sufficient for the system to reach a steady state since the 
local velocity is independent of the site. For simplicity, 
we set to = 7 = Uq = a = 1. 

III. NUMERICAL RESULTS AND DISCUSSION 

Figure 3 (a) shows the average velocity as a function 
of the driving frequency for overdamped case and = 1. 
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FIG. 3: (Color online) Transport of overdamped Brownian 
particles for N — 1. (a) Average velocity v as a function of 
the driving frequency uj; (b) Average velocity v as a function 
of the asymmetry parameter A of the on-site potential for 
uj — 0.1 and uj = 2.0. The other parameters are Ao — 2.0, 
A = 1.0, and D = 0.5. 



(uj = 2.0), one can obtain the opposite velocity, negative 
for A > and positive for A < 0. Moreover, for both 
cases, there is an optimal value of A at which the veloc- 
ity takes its extremal value. When A — > 0, the system is 
absolutely symmetric and directed transport disappears. 
When A — > oo, the asymmetric potential described in 
Eq. (7) reduces to symmetric one with higher barriers, 
U(x) = — ^j-A sin(2a;), resulting in zero velocity. 
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(b) Overdamped case 
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The one-particle stochastic ratchets and stochastic res- 
onance [23|, [24[ has been extensively studied. In the 
adiabatic limit uj — > 0, the ac driving force can be ex- 
pressed by two opposite forces Aq and — Aq. The parti- 
cles get enough time to cross both side from the minima 
of the potential. It is easier for particles to move toward 
the gentler slope side than toward the steeper side, so 
the average velocity is positive. As the frequency to in- 
creases, due to the high frequency, the particles in one 
period have sufficient time to diffuse across the steeper 
side of the well, and therefore it leads to a negative av- 
erage velocity. When the ac driving forces oscillate very 
fast uj — > 00, the particles will experience a time average 

2-7T 

constant force F = J Q " F(t)dt = 0, so the average veloc- 
ity goes to zero. At some intermediate values of uj, the 
average velocity crosses zero and subsequently reverse its 
direction. Figure 3 (b) depicts the average velocity as 
a function of the asymmetry parameter A for different 
frequencies. For the low frequency (uj = 0.1), the aver- 
age velocity is negative for A < 0, zero at A = 0, and 
positive for A > 0. However, for the medium frequency 
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FIG. 4: (Color online) Average velocity v as a function of the 
asymmetry parameter A for TV > 1. (a) Overdamped case: 
uj — 0.1; (b) Overdamped case: uj = 2.0; (c) Underdamped 
case: uj — 0.1 and uj = 2.0. The other parameters are k — 1.0, 
Ao = 2.0, N = 11, and D = 0.5. 



Now we will focus on finding how the interaction be- 
tween the particles affects the directed transport for 
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N > 1. Figure 4 (a) depicts the dependence of the av- 
erage velocity on the asymmetric parameter A for the 
low frequency (w = 0.1). For case A (q = 1), the ac force 
acting on the left end particle, the average velocity is pos- 
itive for A > 0, it is easy for the particles to move from 
the left to the right and the system is a particle conduc- 
tor. However, in the region A < 0, the average velocity is 
zero, the particle can not move from the right to the left 
and the system behaves as a particle insulator. For case 
B (q = N), the ac force acting on the right end particle, 
the transport from the left to the right is prohibited and 
the transport from the right to the left is permitted. The 
average velocity as a function of the asymmetry param- 
eter A for the medium frequency (w = 2.0) is shown in 
Fig. 4 (b). It is found that the unidirectional transport 
will also occur, but the direction of the velocity is oppo- 
site to that in Fig. 4 (a). From Figs. 4 (a) and 4 (b), we 
can conclude that the transport from the free end to the 
ac acting end is prohibited for low frequencies, while the 
transport from the ac acting end to the free end is cut off 
for medium frequencies. Therefore, the system behaves 
as a particle diode. We also found that the similar phe- 
nomena will occur for underdamped case shown in Fig. 
4 (c). 

Note that the similar effects of the particle diode occur 
for the Lennard- Jones potential V(x) — 4fc(^ — 2^), 
where the repulsive force is dominated. However, the ef- 
fects of the particle diode will disappear when the inter- 
action potential is symmetric, for example, the Harmonic 
potential V(x) = ^kx 2 . In order to explain the unipo- 
larity of the transport, we investigate the motion of the 
central particle with zero noises for low frequencies. For 
Morse potential and case A (q — 1), we can find that 
the central particle behaves as a oscillator for A = 1.0. 
In this case, the driving force acting on the first particle 
can be transferred to the other particles. Therefore, di- 
rected transport will occur. However, the central particle 
is static for A = — 1, the central particle can not feel the 
driving force, so the transport is prohibited. For case B 
(q = N), where the driving force is applied at the last 
particle, the central particle is static for A = 1.0 and 
vibrates for A = —1.0. Figure 5 (c) shows the movement 
of the central particle for Harmonic potential and case 
A. We can see that the central particle vibrates for both 
A = 1.0 and A = —1.0, so the system can not behave as 
a particle diode. Therefore, the asymmetry of the trans- 
port originates from the asymmetry of the interaction 
potential. 

Figure 6 (a) describes the average velocity as a func- 
tion of the driving frequency for cases A (q = 1) and B 
(q = N) at A = 1.0. For low frequencies, the transport 
is prohibited for case B. For medium frequencies, the av- 
erage velocity is zero for case A. The average velocity 
for both cases A and B is always zero for the high fre- 
quencies. From Figs. 6 (b) and 6(c), we can see that 
the transport is prohibited for case B at A = 1.0 and 
lu = 0.1. For case A, there exists a value of Aq at which 
the average velocity is maximal. When Aq — )• 0, the ac 
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FIG. 5: (Color online) Time dependent of the position of the 
central particle with zero noises (D = 0.0) for overdamped 
case, (a) Morse potential and q = 1; (b) Morse potential 
and q = N; (c) Harmonic potential and q = 1. The other 
parameters are N = 11, u = 0.1, Aq — 2.0, and k = 1.0. 



driving force disappears, so the average velocity goes to 
zero. When Aq — > oo, the driving force is very large and 
the effect of the potential will disappear, the average ve- 
locity tends to zero. Therefore, there is a peak in the 
curve of v — Aq . Similarly, the curve of v — D shown in 
Fig. 6 (c) for case A is observed to be bell shaped. When 
D — > 0, the particles cannot pass across the barrier and 
there is no directed current. When D — > oo so that the 
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FIG. 7: (Color online) Average velocity v as a function of the 
well depth of the interaction potential k for overdamped case. 
The other parameters are uj = 0.1, Ao = 2.0, A = 1.0, and 
D = 0.5. 



force will disappear, and the velocity also tends to zero. 
Therefore, there exists an optimal value of the well depth 
of the interaction potential at which the average velocity 
takes its maximum. 
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FIG. 6: (Color online) Transport for overdamped case, (a) v 
versus uu at Ao — 2.0 and D — 0.5; (b) v versus Ao at D = 0.5 
and uj — 0.1; (c) v versus D at Ao = 2.0 and ui = 0.1. The 
other parameters are N = 11, A = 1.0, and k = 1.0. 



FIG. 8: (Color online) Dependence of average velocity v on 
the system size N for different cases. For all cases the average 
velocity decrease in the form v oc iV -1 . The other parameter 
are k = 1.0, A = 1.0, and Ao = 2.0. 



noise is very large, the effect of the potential disappears 
and the average velocity tends to zero, also. Therefore, 
one can see that the curves demonstrate nonmonotonic 
behavior. 

The dependence of the average velocity on the well 
depth of the interaction potential is shown in Fig. 7 for 
overdamped case. When k — > 0, the interactions between 
the particles disappears, the central particle can not ex- 
perience the ac driving force, the velocity goes to zero. 
When k — )• oo, the particles are strongly coupled, the 
interacting force is very large, the effect of the driving 



We show the system size dependence of the average 
velocity in Fig. 8 for different cases. It is clearly seen 
that the average velocity decreases monotonically with 
the system size TV for all cases. Remarkably, the average 
velocity decreases in power law v oc iV -1 . Note that the 
average velocity becomes to be independent of the system 
size when the driving force is applied at every particle of 
the chain. 

Since the effect of rectification is very strong, it is nec- 
essary to give a theoretical analysis of the observed phe- 
nomenon by using some approximation. Here, we will 
focus on the small system with N = 3, D = 0.2 and 
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bj = 0.1 for case A. The equations of motion for the two 
particles atfc = a = 7 = l are 

x\ = cos(xi)+^ cos(2xi) — G(x 2 — xi)+£i(t)+A sm(ut), 

(11) 

i 2 = cos(x 2 )+y cos(2x 2 )+G(x2-Xi)-G(x3-X2)+£ 2 (*)> 

(12) 

x 3 = cos(x 3 ) + — cos(2x 3 ) + G(x 3 - x 2 ) + &(*), (13) 

where G(x) = — 2exp(— x)[l — exp(— x)]. 

For case of A = 1.0, when the ac force acts on the left 
end, the particle i = 1 will jump first, and we assume that 
the other particles are oscillating on average about the 
potential minima. Then, Eq. (11) can be approximately 
written 



x'i = cos(xi) + — cos(2xi) — (x2 + x 3 )+£i(t)+Aosin(wt), 

(14) 

the term x*2 + x'3 is a random and does not dominate the 
transport. From this equation we can easily find that the 
particle i = 1 will move towards the right on average. 
The distance between particle i = 1 and i = 2 is small 
and their interaction becomes strong. Then the particle 

1 = 2 can also jump and the ac driving force has been 
transferred from the particle i = 1 to the particle i = 2. 
Therefore, Eq. (12) can approximately become 

x 2 = cos(x 2 ) + cos(2x 2 ) — x" 3 + /isin(u;t) +£ 2 (^), (15) 

where fj, is a constant. From Eq. (15), we can find that 
the particle i = 2 will also move towards the right on av- 
erage. Similarly, the distance between the particle i = 2 
and i = 3 becomes small, and the ac force can be trans- 
ferred to the particle i = 3. The coupling terms prevent 
breaking the chain apart. Therefore, the three particles 
will move together to the right (see Fig. 9 (a)) and the 
particles i = 2 and 3 behave a collective movement. 

However, for A = —1.0, the average velocity based on 
Eq. (14) is negative and the particle i = 1 will move 
to the left on average, the distance between the particle 

2 = 1 and i = 2 becomes large. As we know, it is easy to 
pull the two particles apart for the Morse potential. The 
interaction will disappear quickly and the particles i = 2 
and 3 can not feel the ac driving force from the particle 
i = l. Therefore, Eq. (12) can be approximately written 



x 2 = cos(x 2 ) + y cos(2x 2 ) 



x 3 +&(*). (16) 



From Eq. (16), we can find that the particle i = 2 be- 
haves thermal equilibrium and no directed transport oc- 
curs. Therefore, the particle i = 1 moves to the left on 
average and the particles i = 2 and 3 have zero average 



velocity (see Fig. 9 (b)). The ac driving force breaks the 
chain apart. If the ac driving can be transferred the other 
particles, the all particles can move together to the same 
direction. Due to the asymmetry of the Morse poten- 
tial (it is easy to pull the particles apart), the ac driving 
can not be transferred to the other particle, the ac driv- 
ing particle will escape from the chain and the chain is 
broken apart. 




X -4.0x10 s 




FIG. 9: (Color online) The snapshots of the particles i = 1, 
2, and 3 for different cases, (a) A = 1.0; (b) A = -1.0. The 
other parameters are k = 1, a = 1, Ao = 2.0, ui = 0.1, and 
D = 0.2. 



IV. CONCLUDING REMARKS 

In conclusion, we study the directed transport of inter- 
acting Brownian particles by applying an ac driving force 
at one end of the chain. From the Brownian dynamic 
numerical simulations, we obtain the average velocity of 
the system for both overdamped and underdamped cases. 
One can observe some features not found in the single- 
particle counterparts, for example, the unipolarity of the 
transport. When the low frequency driving force is ap- 
plied at the left end particle (q = 1), the average veloc- 
ity is positive for A > 0, the particles can easily move 
from the left the right and the system works as a par- 
ticle conductor. However, for A < 0, the transport is 
prohibited, the particle can not move from the right to 
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the left and the system behaves as a particle insulator. 
When the low frequency driving force is applied at the 
right end particle (q — N), the transport from the left to 
the right is prohibited and the transport from the right 
to the left is permitted. However, the direction of the 
transport for the medium frequency force is opposite to 
that for the low frequency force. We can call this system 
with the unidirectional transport as the particle diode. 
The unipolarity of the transport is caused by the asym- 
metry of the interaction potential. This asymmetry of 
the transport will disappear in the symmetry interaction 
potential, for example, Harmonic potential. In addition, 



it is also found that there exists an optimal value of the 
well depth k at which the average velocity is maximal. 
The average velocity v decreases monotonically with the 
system size N, v oc TV -1 . 

This work was supported in part by the National Nat- 
ural Science Foundation of China (Grants No. 61072029, 
No. 11004082, and No. 10947166) and the Guang- 
dong Provincial Natural Science Foundation (Grants No. 
01000025 and No. 01005249). Y. -f. H. also acknowl- 
edges the Natural Science Foundation of Hebei Province 
of China (Grant No. A2011201006). 



[1] P. Hanggi and F. Marchesoni, Rev. Mod. Phys. 81, 387 
(2009); P. Reimann, Phys. Rep. 361, 57 (2002). 

[2] F. Julicher A. Adjari, and J. Prost, Rev. Mod. Phys. 69, 
1269 (1997). 

[3] J. Rousselet, L. Salome, A. Adjari, and J. Prost, Nature 
370, 446 (1994); L. P. Faucheux, L. S. Bourdieu, P. D. 
Kaplan, and A. J. Libchaber, Phys. Rev. Lett. 74, 1504 
(1995). 

[4] C. S. Lee, B. Janko, I. Derenyi, and A. L. Barabasi, Na- 
ture 400, 337 (1999); I. Derenyi, C. S. Lee, and A. L. 
Barabasi, Phys. Rev. Lett. 80,1473 (1998); F. March- 
esoni, Phys. Rev. Lett. 77, 2364(1996). 

[5] R. Bartussek, P. Hanggi, and J. G. Kissner, Europhys. 
Lett. 28, 459 (1994); M. O. Magnasco, Phys. Rev. Lett. 
71, 1477 (1993). 

[6] P. Reimann, Phys. Rep. 290, 149 (1997); J. D. Bao and 
Y. Z. Zhuo, Phys. Lett. A 239, 228 (1998);P. Reimann, 
R. Bartussek, R. Haussler, and P. Hanggi, Phys. Lett. A 
215, 26 (1996); B. Q. Ai, L. Q. Wang, and L. G. Liu, 
Chaos, Solitons Fractals 34, 1265 (2007). 

[7] C. R. Doering, W. Horsthemke, and J. Riordan, Phys. 
Rev. Lett. 72, 2984 (1994); R. Bartussek, P. Reimann, 
and P. Hanggi, Phys. Rev. Lett. 76, 1166 (1996). 

[8] B. Q. Ai and L. G. Liu, Phys. Rev. E 74, 051114 (2006); 
B. Q. Ai, Phys. Rev. E 80, 011113 (2009); F. Marchesoni, 
S. Savel'ev, Phys. Rev. E 80, 011120 (2009); B. Q. Ai, J. 
Chem. Phys. 131, 054111 (2009). 

[9] Z. G. Zheng, G. Hu, and B. Hu, Phys. Rev. Lett. 86, 
2273 (2001). 

[10] I. Derenyi and T. Vicsek, Phys. Rev. Lett. 75, 374 (1995); 
F. Julicher and J. Prost, Phys. Rev. Lett. 75, 2618 (1995); 
Z. Csahok, F. Family, and T. Vicsek, Phys. Rev. E 55, 
5179 (1997); A. Igarashi, S. Tsukamoto, and H. Goko, 
Phys. Rev. E 64, 051908 (2001). 

[11] S. Savelev, F. Marchesoni, and F. Nori, Phys. Rev. Lett. 
91, 010601 (2003); S. Savelev, F. Marchesoni, A. Taloni, 
and F. Nori, Phys. Rev. E 74, 021119 (2006); S. Savelev, 
F. Marchesoni, and F. Nori, Phys. Rev. E 71, 011107 



(2005); S. Savelev, F. Marchesoni, and F. Nori, Phys. 

Rev. E 70, 061107 (2004). 
[12] Jose Eduardo de Oliveira Rodrigues and R. Dickman, 

Phys. Rev. E 81, 061108 (2010); W. J. Nie, Q. H. Liao, 

and J. Z. He, Phys. Rev. E 82, 041130 (2010). 
[13] M. Kruger and M. Fuchs, Phys. Rev. E 81, 011408 (2010); 

M. Evstigneev, S. von Gehlen, and P. Reimann, Phys. 

Rev. E 79, 011116 (2009); P. Reimann, R. Kawai, C. 

Van Den Broeck, and P. Hanggi, Europhys. Lett. 45, 545 

(1999). 

[14] R. M. da Silva, C. C. de Souza Silva, and S. Coutinho, 
Phys. Rev. E 78, 061131 (2008); A. D. Chepelianskii, M. 
V. Entin, L. I. Magarill, and D. L. Shepelyansky, Phys. 
Rev. E 78, 041127 (2008). 

[15] G. Szamel, J. Chem. Phys. 127, 084515 (2007). 

[16] A. E. Filippov, J. Klafter, and M. Urbakh, Phys. Rev. 
Lett. 92, 135503 (2004); F. Family, H. G. E. Hentschel, 
and Y. Braiman, J. Phys. Chem. B 104, 3984 (2000). 

[17] E. Heinsalu, M. Patriarca, and F. Marchesoni, Phys. Rev. 
E 77, 021129 (2008); X. R. Qin, B. S. Swartzentruber, 
and M. G. Lagally, Phys. Rev. Lett. 85, 3660 (2000). 

[18] C. Lutz, M. Kollmann, and C. Bechinger, Phys. Rev. 
Lett. 93, 026001 (2004). 

[19] K. Luo, T. Ala-Nissila, S.-C. Ying, and A. Bhattacharya, 
Phys. Rev. Lett. 100, 058101 (2008) 

[20] L. L. Bonilla, Phys. Rev. B 35, 3637 (1987). 

[21] H. S. J. van der Zant, T. P. Orlando, S. Watanabe, and 
S. H. Strogatz, Phys. Rev. Lett. 74, 174 (1995). 

[22] R. L. Honeycutt, Phys. Rev. A 45, 600 (1992). 

[23] V. N. Smelyanskiy, M. I. Dykman and B. Golding, 
Phys. Rev. Lett. 82, 3193 (1999); M.I. Dykman, H. Ra- 
bitz, V. N. Smelyanskiy, B.E. Vugmeister, Phys. Rev. 
Lett. 79, 1178 (1997); V. N. Smelyanskiy, M. I. Dykman, 
H. Rabitz, and B. E. Vugmeister, Phys. Rev. Lett. 79, 
3113(1997). 

[24] D. G. Luchinsky, M. J. Greenall, P. V. E. McClintock, 
Phys. Lett. A 273, 316 (2000). 



